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Long Range Tensor Correlations in 
Charge and Parity Projected Fermionic Molecular Dynamics 

Sonia Bacca*, Hans Feldmeier and Thomas Neff 
Gesellschaft fiir Schwerionenforschung, Planckstr, 1, 64291 Darmstadt, Germany 

Within the framework of Fermionic Molecular Dynamics a method is developed to better account 
for long range tensor correlations in nuclei when working with a single Slater determinant. Single- 
particle states with mixed isospin and broken parity build up an intrinsic Slater determinant which 
is then charge and parity projected. By minimizing the energy of this many-body state with respect 
to the parameters of the single-particle states and projecting afterwards on angular momentum 
ground state energies are obtained that are systematically lower than corresponding Hartree-Fock 
results. The realistic Argonne V18 potential is used and short range correlations are treated with 
the Unitary Correlation Operator Method. Comparison with exact few-body calculations shows 
that in *He about one fifth of the correlation energy due to long-range correlations are accounted 
for. These correlations which extend over the whole nucleus are visualized with the isospin and 
spin-isospin density of the intrinsic state. The divergence of the spin-isospin density, the source for 
pion fields, turns out to be of dipole nature. 

PACS numbers: 21.60.-n, 21.60.Do, 21.10.Hw, 21.30.Fe 
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I. INTRODUCTION 

Despite the undeniable successes of mean field mod- 
els a more microscopic view on the structure of atomic 
nuclei reveals correlations among the nuclcons of various 
kinds. These arc induced by realistic nucleon-nuclcon in- 
teractions and cannot be represented by a single Slater 
determinant in a mean field picture. 

Impressive progress has been made in the derivation of 
the nuclear interaction, one of the latest achievements 
being the nucleon-nucleon (NN) forces and consistent 
many-nucleon forces developed with chiral perturbation 
theory [l|, [2]. Furthermore, a number of other nuclear 
potentials exist, like the scmi-phenomenological Argonne 
V18 potential (AV18) Q, which reproduce the experi- 
mental NN phase shifts with high precision. Common to 
the realistic NN forces mentioned before are two general 
features: (i) they are strongly repulsive at very small dis- 
tances, preventing nucleons to stay close together, which 
induces short-range repulsive radial correlations among 
them; (ii) they contain a tensor force component, such 
that the nucleon pair feels a force which depends on 
their spin orientation with respect to the relative dis- 
tance. This induces further short and long range tensor 
correlations among nucleons. 

It is well known that a Slater determinant (antisym- 
metrized product state) cannot represent such corre- 
lations. Models that use a Slater determinant basis 
need very large many-body Hilbert spaces to represent 
these correlations. Therefore in the no-core shell model 
(NSCM) U the unitary Lec-Suzuki transformation is in- 
corporated to treat the short range part of the Hamil- 
tonian that scatters to very high lying oscillator shells 
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and thus helps to improve convergence and to reduce the 
dimensions of the Hilbert space. 

Another possible solution to overcome the problems 
caused by the short-range correlations is the Unitary Cor- 
relation Operator Method (UCOM) 0, H, [3|- A unitary 
correlator C = C^Cr is devised to imprint explicitly the 
short range tensor correlations, by means of Co, and 
radial correlations, by C^, into the uncorrclatcd Slater 
determinant basis. 

However, tensor correlations are not only of short range 
like the radial repulsive correlations, but contain parts 
of long-range nature originating from the one pion ex- 
change, exemplified e.g. by the large extension of the D- 
wavc component of the dcuteron wave function. There- 
fore, even when using the UCOM approach long-range 
correlations need to be incorporated by configuration 
mixing in the Hilbert space. In case of many-body meth- 
ods based on an expansion of the states in terms of a 
complete set of basis states one is in principle able to 
account for all kind of correlations, provided that con- 
vergence in the expansion is reached. But when working 
with a low-momentum basis, like with just one or few 
Slater determinant states, the incapability of describing 
long-range correlations in nuclei still constitutes a big 
limitation. 

In this work, we would like to address the general 
problem of how to account for long-range tensor cor- 
relations within a restricted Slater determinant basis. 
We propose a method, which follows ideas presented by 
Ikeda, Sugimoto and Toki published in a scries of papers 
[3, [lO, llll, [l3, 13 to allow proton and neutron states 
to mix at the single-particle level and then project the 
many-body state on good charge and parity. We in- 
troduce additional variational degrees of freedom in the 
Fermionic Molecular Dynamics (FMD) framework which 
allow for the possibility that the isospin of a nucleon can 
point in any direction in isospin space and not just in the 



proton or neutron direction. 

The paper is organized as follows. In Sec. [TT] the theo- 
retical background is set: an overview of the FMD and of 
its extension and the used effective interaction are pre- 
sented. In Sec, mil results are shown and conclusions are 
finally drawn in Sec. IIVI 



II. THEORETICAL OVERVIEW 

A. Fermionic Molecular Dynamics 

In the following we will briefly introduce the Fermionic 
Molecular Dynamics (FMD) approach and the gener- 
alization of the FMD wave functions to mixed charge 
states. In FMD [13, M, M, [13, M, [II the many-body 
Hilbert space is spanned by non-orthogonal many-body 
basis states which are given by antisymmetrized products 
of single particle states, as 



Q) = T^(ki)®-"® k^)) 



(1) 



A is the antisymmetrization operator and q^ ) de- 
note the single-particle states which are linear combina- 
tions of Gaussian wave packets localized in phase-space, 
with variable spin orientation x ) ^-nd the generalized 
isospinor \(,), 



( -^ k ) = X] '^* '^^P 1 



2a, 



®\x^)®\£.i) ■ (2) 



Each Gaussian wave packet is parameterized in terms 
of a complex vector bi , indicating the mean position and 
momentum and a complex width parameter a^, which can 
be different for each Gaussian, in contrast to the AMD 
approach [20] where the widths are real and common for 
all nucleons. 

The spinor is parameterized via two complex compo- 
nents for spin-up and spin-down 
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allowing for all orientations of the spin. Analogously the 
generalized isospin part | ^i ) describes a linear superpo- 
sition of proton p j and neutron n ) 



|cO = ^np)+er>) 



(4) 



so that a nuclcon can adopt any "direction" in isospin 
space. In former applications of FMD | S, ) was either a 
proton p ) or a neutron state Ui ) . 

This generalization introduces charge mixing in the 
single-particle and many-body space. In a description, 
where the degrees of freedom are only nucleons and the 
charged mesons do not appear explicitly, the charge car- 
ried by the nucleons is a sharp quantum number. The 
many-body state Q ) however breaks the symmetry with 



respect to isospin rotations around the 3-axis. Therefore 
one has to project on the desired charge number corre- 
sponding to the eigenvalue Mt of the third component 
of the total isospin 
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This is achieved by the charge projection operator [2^ 
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Thus, the charge projected state, can be written as a su- 
perposition of Slater determinants | Q^"^ ) obtained by 
rotating the single determinant \Q) in isospin space 
about an angle 2Trn/A: 
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The following simple example shows that the projec- 
tion of a two-body product state with mixed charges re- 
sults in correlated states with good isospin. From the 
product states of two nucleons 



(I 



p)± jn)") (K) (\p) + \n)^ = 

\p) ®\p) ± |n)(g)|n) 

+ (\p)(g)\n)±\n)(g)\p)'\ (8) 

one can project out all 4 components with isospin T =1 
and including the Mt = two-body states 

\T = 1,Mt = Q) = -^(\^p) ®\n) + \n) ®\p)^ (9) 

|T = 0,AfT = 0) = —{\^p) ®\n) -\n) ®\p)^ (10) 

which are correlated and cannot be written as a product 
of two single-particle states. 

As the lightest charged meson, the pion, is of pseu- 
doscalar nature with a negative intrinsic parity it may 
carry besides charge also parity from one nuclcon to the 
other. Therefore we also allow for parity breaking in the 
FMD state, which simply means that the parameters hi 
are not the same for all nucleons, and restore it by projec- 
tion on good parity 7r = ±l with the projection operator 
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(11) 



where 11 is the parity operator. 

In order to sec in how far Slater determinants, that 
break charge and parity, can represent long range cor- 
relations induced by the exchange of pions variational 



calculations arc performed. Wc minimize the energy 

( Q; t:, Mt | iJ | Q; tt, A/t ) {Q\H P^^^P^ \ Q ) 

(Q;7r,MT|g;7r,AfT) " (Q|P^^P-|Q) 

(12) 
of the charge and parity projected FMD state 

\Q;tt,Mt)=P^'^P^\Q) (13) 

with respect to all single-particle parameters contained 
in the single Slater determinant Q ) • 

As the intrinsic Hamiltonian H commutes with 11 and 
T3 one has to project only the ket-state so that the num- 
ber of terms in the energy (|12p is reduced. 

The correlated many-body state I Q;tt,Mt) that re- 
sults from the minimization of the energy ([T^ in general 
breaks rotational and translational symmetry. Therefore 
we project after the variation on good angular momen- 
tum and center of mass momentum zero |18l . Il9j : 



l,r,MK,AlT) 



PilKPcM\Q;7T,MT) 



(14) 



Throughout this paper we consider only energies of J^ = 
states so that K = M = 0. All energies are calculated 
as expectation values of the intrinsic Hamiltonian H = 

T — Tcm + y where the center of mass kinetic energy has 
been subtracted: 



E{0+) 
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It would be of course desirable to first restore all sym- 
metries of the Hamiltonian and then do the variation 
with respect to the single-particle parameters contained 
in \Q)- However, angular momentum projection re- 
quires the superposition of a few hundred rotated states, 
which is numerically very costly. Therefore we do a vari- 
ation after charge and parity projection and project on 
angular momentum and CM-momcntum zero after vari- 
ation. 



B. Unitary Correlation Operator Method 

Realistic nucleon nucleon interations, like the chiral 
forces or the Argonne V18 potential, that reproduce the 
phase shifts and the deuteron properties, induce strong 
short range correlations. The repulsive core prevents two 
nuclcons from getting too close. In the T = channel the 
tensor force induces correlations by aligning the spins 
of two interacting nucleons along their distance vector. 
Both correlations cannot be represented by Slater deter- 
minants. Even though the Slater determinants of a shell 
model basis form a complete set, diagonalization of a re- 
alistic Hamiltonian in the 4-body space shows that con- 
vergence cannot be reached in any tractable shell model 
space [8|. Thus many-body Slater determinants are an 



inadequate representation for these short ranged corre- 
lations. In fact the no-core shell model employs the 
Lee-Suzuki transformation of the Hamiltonian which im- 
proves the convergence dramatically. 

Here we will use the Unitary Correlation Operator 
Method (UCOM) to take care of the short range cor- 
relations. The reason is that an effective interaction to 
be used in the FMD Hilbert space has to be represented 
in an operator form, rather than in terms of matrix ele- 
ments as is the case for a Lee-Suzuki transformation or a 
G-matrix. UCOM provides both representations so that 
we can compare with no-core shell model results. 

The concept of UCOM 0, [y, Q consists in treating ex- 
plicitly the strong short range correlations induced by the 
hard core and the short range part of the tensor force by a 
state independent unitary transformation, while the long 
range correlations have to be represented by the many- 
body states spanning the Hilbert space. When the cor- 
relation operator is applied to an initial Hamiltonian, a 
phase-shift equivalent correlated interaction is obtained. 
The correlated Hamiltonian is defined via a similarity 
transformation 
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where Cr and Cn are the unitary central and tensor cor- 
relation operators, respectively. Cr shifts pairs of nu- 
cleons radially away from each other whenever their dis- 
tance is so small that they would be inside the repulsive 
core. Cn aligns nucleon pairs with T = and S^l along 
the direction of their total spin so that the typical tensor 
correlations known from the deuteron are imprinted into 
the many-body state. Different from the radial correla- 
tions induced by the short ranged repulsion, the tensor 
interaction mediated by the exchange of pions, the light- 
est of all mesons, induces long range correlations. In 
order to keep the effect of the induced many-body inter- 
actions small one restricts the range of the action of the 
tensor correlation operator Cq and represents the long 
range part of the tensor correlations with the many-body 
state (see Ref. Q for details). 

The two-body part of the correlated Hamiltonian (fT6|) 
is used as an effective intrinsic Hamiltonian 



H = T -Tr 
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UCOM 



(17) 



that is applicable in low momentum Hilbert spaces. The 
effective interaction Vucom is phase shift equivalent to 
the initial realistic interaction. Moreover different real- 
istic potentials, like Bonn A, Nijmegen or Argonne V18, 
lead to practically the same Vucom- In FMD the effec- 
tive potential needs to be given in operator representa- 
tion and not as momentum space matrix elements for the 
different partial waves. Vucom can easily be obtained in 
operator form if the initial interaction is given in terms of 
operators. Therefore we use in this publication the Ar- 
gonne V18 potential (AVIS) as initial potential because 
it is already in operator form. 



For a better understanding of the results let us first 
consider in somewhat more detail how the tensor corre- 
lator Cfi acts. Any operator that depends only on the 
distance r between two nucleons is invariant under ac- 
tions of the tensor correlator 



Ci,'V{r)C, 



nz) 
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because the tensorial generator 
3 



gn = "^{r)-{{pn-ai){r-G2) + (^•cti)(po-CT2) ) (19) 

commutes with r ~ \r\. The reason is that the so called 
orbital momentum pn, defined as the component of the 
relative momentum perpendicular to r, 



on mean-field theories like Hartree-Fock or simple shell- 
model pictures. The reason is that those calculations can 
be performed without the tensor force using only central 
and spin-orbit interactions that are fitted to energies and 
other properties of the many-body system. The success 
of these models does not imply that in finite nuclei or 
nuclear matter the tensor force is not important. By 
fitting the parameters of the interaction one has effec- 
tively moved the correlation energy of the tensor into 
central and spin-orbit interactions. However, the result- 
ing effective interactions are not phase-shift equivalent 
any longer, Vucom is. Furthermore, the uncorrelated 
many-body Slater determinants do not possess the ten- 
sor correlations anymore that are existent in reality and 
can be observed for example in the high momentum part 
of the momentum distribution of the nucleons [dl . 



Pn =P 



(20) 



III. RESULTS 



commutes with r. The correlation function ^(r) defines 
the strength of the transformation as a function of the 
distance between an 5*=! nucleon pair. 

In the following we shall use different d{r) (displayed 
in Fig. [T]) with varying range characterized by the range 
parameter I^ ~ j dr r'^'d{r). 

When applying the tensor correlator to the tensor in- 
teraction 

VT(r)S,2 - VT(r) ( A(r.ai)(r.52) - (£i-£2)) (21) 



In the following subsections we study in how far long 
range correlations originating from the tensor interaction 
can be repesented by a charge and parity projected Slater 
determinant. For that the realistic AV18 interaction is 
used and the results are compared to exact solutions for 
the 4-body system "^He. 

In subsection IIII CI we use a phenomenological inter- 
action that is not based on a realistic nucleon nucleon 
force to illustrate that fitting interaction parameters to 
a specific nucleus may be misleading when drawing con- 
clusions about the strength of correlations. 



it is unitarily "rotated" into a reduced tensor Vr, a cen- 
tral Vc, and a spin-orbit Vls component, plus some other 
terms which are negligibly small. 



Co' VT{r)Su Cn = Vt + Vc + Vls + 



(22) 
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In Eq. (P^ one readily sees that the induced central 
force Vc increases if the correlation function t9(r) in- 
creases. Thus whith larger range of the tensor correlator 
more strength goes to the central and spin-orbit inter- 
action. One should keep in mind that the correlation 
function i?(r) is different for isospin T = and T = I, 
see Ref. [y| and that there are also contributions from the 
correlated kinetic energy C^^^ T C^, the spin orbit force 

Co 7-s Cq and possibly other terms in the interaction. 



Eq. ((22)) also gives finally an explanation why in most 
textbooks after the discussion of the deuteron with its 
quadrupole moment and the mandatory tensor force the 
tensor interaction dissappears in the following chapters 



A. Correlation energies 

In this section we investigate how the induced corre- 
lations are accounted for by breaking charge and parity 
in each single-particle state as a function of the strength 
of the tensor component of the nucleon-nucleon interac- 
tion. In order to get a fair judgement we use the realistic 
effective interaction Vucom which is derived from the 
AV18 potential. The strength of its tensor component 
depends on the range I^ of the tensor correlator. This 
range corresponds to the cutoff in momentum space in 
the Mow fc-approach 0, [23]. For all values of I^ the effec- 
tive interaction Vucom is phase shift equivalent to the 
original realistic interaction. It should be noted that for 
different I§ all parts of the Vucom potential change ac- 
cordingly under the similarity transformation (J16p and 
not just the tensor part (see also Eq. (HH)). 

We calculate the ground state energy of ^He in three 
different ways using the same effective Hamiltonian p7)) 
for different I^. 

First we minimize the energy (Q|-ff|(3)/(Q|Q) 
with respect to all FMD parameters contained in Q with- 
out any projection. The FMD single-particle states I q^ ) 
contain two Gaussians, see Eq. ^. The results are la- 
belled HF as this variation corresponds to a Hartree-Fock 
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FIG. 1: (Color online) Upper panel: tensor correlation func- 
tions i9(r) with increasing range I^ (arrow). Lower panel: 
ground state energies of ''He calculated with the HF and CPP 
states compared to exact results. The correlation functions 
shown in the upper panel are used. 



minimization. Second we minimize the energy (|12p of 
the charge and parity projected FMD Slater determinant 
(fT3ll . labelled CPP. 

After that an angular momentum and CM projection 
is applied to the HF and CPP state, sec Eq. dHD. The HF 
state has already J'^ ~ 0+ and Mt ~ but the CPP state 
has no good angular momentum. The resulting energies 
P^ are compared to the exact results of no-core shell 
model and hyperspherical harmonics calculations. 

In order to study the effects of different tensor 
strengths in the Hamiltonian we vary the range of the 
tensor correlator in the 5* = 1, T = channel, which ac- 
cording to Eq. (|22p moves strength from the tensor to the 
central part. This has the advantage that we are always 
using phase-shift equivalent potentials. In the upper part 
of Fig. [T] we display the correlation functions employed. 
They differ in range which we quantify by the range pa- 
rameter I^. 

In the lower part of Fig. [1] we present the binding en- 



ergies of ^Hc obtained with the Vucom potential for dif- 
ferent I^ in the three variational Hilbert spaces as ex- 
plained before. The range Ig = 0.09 fm"^ corresponds 
to the standard choice where Vucom reproduces within 
0.5 MeV both the '^He and '^He binding energies. This 
interaction has been used for many apphcations, ranging 
from few- to many-body problems [i,t2ll,[2l, [11,12^. The 
range /^ — 0.09 fm is optimal in the sense that the net 
contribution from 3-body forces is minimized [8| . 

The exact results in Fig. [T] show that ranges smaller 
than 7^=0.09 fm"^ lead to undcrbinding with respect to 
the experimental energy of -28.3 MeV (like for most 
realistic potentials). On the other hand tensor corre- 
lators with ranges larger than 7^=0.09 fm"^ induce 3- 
and 4-body interactions in the correlated Hamiltonian 
H, Eq. (|16p . which are on average repulsive but not in- 
cluded here. A similar kind of overbinding was obtained 

1.6 fm~ 
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The first result to be observed in Fig. [T] is the decreas- 
ing difference between the HF energy and the exact result 
when the range /^ of the tensor correlator is enlarged. 
This can be understood if one keeps in mind that the 
HF state of '*Hc is a pure (Os)^ configuration so that the 
tensor interaction cannot contribute. All binding comes 
from the central part of the interaction. According to 
Eq. (P^ increasing the correlation strength implies that 
the induced central part of the correlated tensor inter- 
action (1221) in the 5=1 r = channel increases and 
binding energy is gained, while the long range part of 
the tensor force, not seen by the HF state, is reduced. 
The difference between the HF and the exact energy is 
the correlation energy which is due to the interaction in- 
duced correlations present in the exact many-body state 
but absent in the HF Slater determinant. 

The exact results obtained in an "unrestricted" Hilbert 
space show much less variation as function of I^ be- 
cause there the long range correlations are represented 
in the many-body state even for small I^. If we would 
include the induced 3- and 4-body potentials the exact 
energy would not depend on I^ because the transforma- 
tion would then be unitary even in the many-body space. 

An important finding is that, although the variational 
manifold {Q} includes the new isospin mixing degrees of 
freedom, the Hartree-Fock (HF) minimum does not make 
use of them. All ^f and ^f are either or 1, so that the 
Hartrec-Fock state is already an eigenstate of charge. It 
turns out to be also an eigenstate of angular momentum 
and parity (J'^ = 0+). 

This changes when the variation is performed after 
charge and parity projection. Now the ^f and ^" pa- 
rameter assume values different from or 1. As Fig. [1] 
shows, the charge and parity projected state (CPP) can 
represent part of the long range correlations leading to a 
correlation energy of about 20% the full correlation en- 
ergy. 

In Table |T] the expectation values of the total Hamil- 
tonian, the kinetic energy, the interaction energy and 
the tensor part of the correlated interaction Vucom are 



TABLE I: Expectation values of different terms of the Hamil- 
tonian for the ^He J'^ = 0+ ground state. The variation of 
single-particle states is performed for Hartree-Fock (HP) , Par- 
ity Projected (PP) and Charge and Parity Projected (CPP) 
intrinsic Slater determinants. Vucom with tensor correlator 
range /^ = 0.09 fm^ is used. Numerical values are in MeV. 



HF[0+] 
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-16.42 


-17.51 


-18.61 


il) 


50.25 


51.56 


57.93 


(Y) 


-66.68 


-69.07 


-76.54 


(Yt) 


0.00 


-1.35 


-4.59 



listed. The HF state has zero tensor energy and the 
smallest kinetic energy, while the CPP state contains 
correlations that yield about —2.2 MeV more binding. 
These consist of 7.7 MeV kinetic energy counterbalanced 
by —9.9 MeV additional potential energy. Out of that the 
tensor contributes —4.6 MeV. The rest originates from 
other parts. The rather large increase in kinetic energy 
is partly due to the fact that in the CPP minimum the 
two Gaussians in the single-particle states are displaced 
from each other, thus break parity and include higher an- 
gular momenta. Therefore the spatial and spin parts are 
different in the CPP minimum and the HF minimum. 

It is interesting to note that a variation after parity 
projection (PP) already takes into account part of the 
correlations as can be seen from the center column of 
Table HI One should also keep in mind that with respect 
to an uncorrclatcd Slater determinant a by far larger 
amount of correlation energy resides in the short range 
tensor correlations which are treated here by the Uni- 
tary Correlation Operator Method (UCOM). These short 
range correlations are so strong that they even cannot be 
properly accounted for in a very large scale shell model 
basis, not to mention a single Slater determinant. 

The effects are not as pronounced if one restricts the 
single-particle states to one Gaussian. Adding the second 
Gaussian is essential in describing tensor correlations, as 
it allows to add higher angular momenta components to 
the single-particle wave functions. 

One should pay attention to the fact that additional 
correlations in the many-body state that are induced by 
the two-body interaction give an attractive contribution 
from the potential but at the same time the kinetic energy 
rises. For example tensor correlations imply admixtures 
of higher angular momenta and thus more kinetic energy. 
Hence the correlation energy is the result of a subtle in- 
terplay between enhanced kinetic energy and increased 
attractive potential energy. 

The ground state energies labelled "Exact" are calcu- 
lated within the no-core shell model (NCSM) using the 
translationally invariant harmonic oscillator formulation 
of Petr Navratil [20] and the hyperspherical harmonics 
(HH) approach developed by Nir Barnea [30, |3l|. The 
latter method has been applied to the calculation of ex- 
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FIG. 2: (Color online) Convergence patterns for the NCSM 
calculations with the UCOM potential for J^ = 0.06 fm'^: the 
ground state energy of He is plotted as a function of the 
oscillator parameter hO, for different model space sizes (upper 
panel) . Comparison of the NCSM and HH convergence as a 
function of A'^max and -R'max, respectively, for fixed fiQ (lower 



panel) . 



act clectroweak reactions on light nuclei mainly with lo- 
cal interactions [H, [H, [13, [S S 113, HI . Recently it 
has been extended to use interactions represented in har- 
monic oscillator (HO) states (see e.g. Refs. [20, [39|). 
The HH approach is equivalent to the NCSM, but makes 
use of HH functions instead of HO eigenstates. The label 
"Exact" is meant in the sense that one systematically en- 
larges the Hilbert space and controls the convergence to 
approach the final many-body state which should contain 
all kind of correlations induced by the interaction. 

Convergence patterns for the NCSM calculations for a 
specific value of I^ = 0.06 fm'^ are shown in Fig. [51 For a 
given size of the model space, characterized by the maxi- 
mum oscillator quantum number A^max, the ground state 
energy is plotted as a function of the oscillator param- 
eter hfl. Convergence, resulting in a flat energy curve 
over a significant range of oscillator parameters, can be 



TABLE II: Exact ground state energies (in MeV) of ''He for different ranges of the tensor correlator in UCOM, calculated with 
the NCSM and HH approaches. 



h [fm=*] 


0.06 


0.09 


0.12 


0.15 


0.20 


0.30 


NCSM 
HH 


-26.80(3) 
-26.84(6) 


-27.80(3) 
-27.84(7) 


-28.62(9) 
-28.62(7) 


-29.26(13) 

-29.25(8) 


-29.86(10) 
-29.92(10) 


-30.70(6) 
-30.86(10) 



obtained already with A^max = 18. The energy gained 
compared with the results obtained with more moderate 
model spaces can be attributed to the residual long range 
correlations not described via the UCOM. 

In Fig. [2] we also compare the convergence of the HH 
and NCSM calculations for two fixed values of hfl as a 
function of the A',„ax parameter, which is for the HH the 
analo g to iV^ax for the NCSM. As already pointed out in 
Refs. [2a] and [40] the convergence of HH is superior to 
the one of the NCSM, since no hil dependence is observed 
even for small Hilbert spaces. Therefore we employ an 
exponential fit of the NCSM energies as a function of 
TVniax (for fixed hU) to extrapolate to infinite dimensions. 
The further energy gain is of the order of 50 — 100 keV. 
The energies obtained for different values of the tensor 
correlator volume Lg used in the UCOM are shown in 
Table [III The NCSM and HH results nicely agree with 
each other within the error bars. 

The small difference of about 0.5 MeV we find in our 
exact calculations for I^ = 0.09 fm^ with respect to the 
values previously published in Refs. [g] with the NSCM 
and [26| with HH for the same I^ is related to some minor 
differences of the potential used. They originate from 
the fact that in this paper we adopt a Vucom potential 
in operator representation which is suited for the FMD 
code. This implies a truncation of the Baker-Campbell- 
Hausdorff expansion (for details see Ref. [3]), not needed 
if one correlates directly the two-body states. 

Here, we would like to stress again that the exact ap- 
proaches make use of thousands of basis states to reach 
convergence in energy, whereas with our improved FMD 
wave function we vary the parameters of a single Slater 
determinant projected on charge and parity. 



B. Pseudo-scalar iso- vector correlations 

A single Slater determinant can usually not represent 
correlations except those induced by the Pauli principle. 
However, if the Slater determinant represents an intrinsic 
state the situation is different. The physical state that 
has the same symmetries as the Hamiltonian is obtained 
by means of projection on angular momentum, parity 
and, as in our case, on charge. It is in general not a 
single Slater determinant but a superposition of many, 
see Eqs. ([T^[H)) . with the restriction that all these Slater 
determinants are generated from a single one by means 
of rotations, parity inversion or rotation in isospin space. 

A well known example for long range correlations in 



nuclei are intrinsically deformed nuclei where out of one 
deformed intrinsic state one projects a whole rotational 
band. 

The intrinsic state for ^He obtained in this paper by 
variation after parity and charge projection shows a pro- 
nounced long range correlation in the spin isospin degrees 
of freedom. To elaborate on that let us consider the ex- 
change of the pseudo-scalar iso- vector pion. 

In all realistic nucleon-nucleon interactions the one- 
pion exchange is responsible for the longe range tail of 
the potential. This part is not affected by the unitary 
correlator Co which is of short range. The induced long 
range correlations should therefore be represented by the 
many- body state. 

The vertex describing the interaction of a nucleon field 
A^(x) with a pion field <f>^ (x) has in pseudo- vector cou- 
pling the following form 



Cn-k{^) 



9-n 

2M 



^iV(x)757^T«7V(x)a^$«(x) , 



i=l 



(23) 

where r'^'^, i = 1, 2, 3 denote the Pauli matrices in isospin 
space. 



r(l): 



1 

1 



r(2) 




In the stationary case the three components of the pion 
field $^ (2) satisfy the time-independent Klein-Gordon 
equation 

(-V^+ m^)$W(.5) = ^V.§^'^5), (25) 

where the source term is the divergence of the nuclear 
isospin current density 



S^"\x)^\N[x)-i''^T('^N{x) 



(26) 



'.ii), 



After a non-relativistic reduction to leading order S (x) 
becomes the one-body spin-isospin density of the nuclear 
many-body system 

FH5) = \Y,^%l-r{k)) a{k) r(^Hk) . (27) 

fc=i 



The relation with the physical pion fields are 

n-{x) = ^{^i'\x) - t<^i^\x)) . (28) 

^°(5) = $?^(x) 

The (3)-component S (x) is the difference between 
the proton and neutron spin density at position x and 

— —(3) 

\/-S (x) is the pseudo-scalar iso- vector source density 

for the 7r° field, while V-S (5) and V-S (x) arc the 
sources for tt* fields. 

In Fig. [3] various intrinsic densities are displayed as a 
function of x and y at the plane z = 0. All densities are 
calculated with the parity projected intrinsic state 



Q; + ) 



p 



7r=+l| 



') 



(29) 



that was obtained by minimizing the energy after parity 
and charge projection (CPP), cf. Eq. P^ . 

Before we come to the spin-isospin density that is 
related to the pion fields we show in Fig. [3] the pro- 
ton and neutron densities for point-like nuclcons. Both 
are equal up to very small deviations caused by the 
Coulomb interaction. This means that the isospin den- 

f 31 

sity pr {x) = Pp{x) — pn{x) in the isospin (3)-direction 
is zero. However, the isospin density 



pW(5) = 



{Q-, + \Eti^H^~r{k))r('Hk)\Q; + ) 



(Q; + |Q; + ) 



(30) 



in the (l)-direction orthogonal to the proton/neutron (3)- 
direction assumes a non-zero value. As seen in Fig.[3]this 
density is of quadrupole type. The intrinsic state Q; + ) 
is apparently not an eigenstate of T^^> and its new de- 
grees of freedom, S,^, £,""■, that allow to mix protons and 
neutrons, are responsible for this density. The charge and 
parity projected state \Q;tt = +1,Mt = ) has (like the 
Hartree-Fock state) again a vanishing isospin-(l) density. 
The lower panels in Fig. [3] display the spin-isospin den- 
sity 



S {x) = 



+ ||(i)(2)|0; + ) 



(31) 



=,(1), 



and its divergence V ■ S {x) which is the source density 

for the $Jr (x) pion field. S (x) represents a pseudo- 
vector iso-scalar field ((l)-component) with a pronounced 
dipole shape. The divergence is the according pseudo- 
scalar iso-vector source density. One should note that 
the structure of the intrinsic state extends over the whole 
nucleus and is hence of long range. 

All other spin-isospin densities are two orders of mag- 
nitude smaller which means zero within numerical un- 
certainty and hence not displayed. One should keep in 



mind that the intrinsic state | Q;-l-) can be rotated in 
isospin space around the (3)-axis resulting only in an 
overall phase of the ground state, because the intrinsic 
state is projected on good charge number by summing up 
rotations around the (3)-axis, cf. Eq. ([6]). Likewise one 
can rotate the intrinsic state in coordinate space with- 
out affecting the angular momentum projected 0"*" ground 
state, so that the dipole in y-direction could also point 
in any other direction. 

Let us try to explain the physical meaning of the 
non-zero intrinsic pseudo-scalar isovector source density 

V-S (x) with help of an analogy to the Coulomb in- 
teraction. Consider a positronium, negativly charged 
electron plus positivly charged positron, in their atomic 
ground state. This state has angular momentum zero and 
the probability to find a positron at some position equals 
that of the electron. Hence the mean value or expectation 
value of the charge density pe is zero and consequently 
there is no Coulomb field $e which of course must not be 
interpreted that there is no Coulomb attraction. But if 
we take a "snap shot" we find the positron and electron 
on opposite sides (perfect correlation) forming a dipole 
with non-zero charge density and non-zero Coulomb field 
by which they attract each other. After projecting this 
dipole on a 0"'' state we get the true positronium ground 
state. 

Looking at Eq. (QH) and replacing 



M 



Airpe 




one recovers the well known equation for a Coulomb 
field created by a charge density. The analogy is obvi- 
ous, the intrinsic state is the "snap shot" where one sees 
the dipole like source density as a one-body mean field 
(Fig. [3]) . If one wants to see that correlation in the spher- 
ical quantum state with good parity one would have to 
resort to two-body information or correlation functions, 
cf. Fig. (2) inRef. 0. 



C. Phenomenological interaction 

In Refs. [§, M, HH, [H, \M Ikeda, Sugimoto and Toki 
propose the idea to mix proton and neutron wave func- 
tion at the single-particle level and then perform charge 
and parity projection of the many-body state. They use 
for the 4 single-particle states for ^He (t/ = 1,2; m ~ 
±1/2) 

(5|j/,im)= J2 (l''^,iraXx)[Y\9,(b)\^)]l®\mt) . 

i = 0,l 

Here, the spatial part of the wave function contains s- 
and p-wave components by construction, and the radial 
part 4>v,imt is expanded in terms of Gaussian functions. 
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FIG. 3: (Color online) Densities calculated with parity projected intrinsic state | Q; + ) obtained for /^ — 0.09 fm^. 

Upper panel; proton density Pp(iE) and neutron density /9„(3) (in fm""^); lower panel: isospin density pr [x) (in fm~^) 

spin-isospin density S (S); divergence of spin-isospin density V ■ S (x) (in fm~*). Contour lines in decades. 



A phenomenological interaction that is composed of a 
central Volkov potential [4l| and the tensor plus spin- 
orbit GSRS force [43] is used. It is argued that the mixing 
of parity and isospin is able to account for appropriate 
tensor correlations in the a-particle when the triplet-even 
component of the central part is reduced by 0.81 and 
the TiTj part of the tensor force enhanced by a factor 
1.5. For this special interaction, which we will refer to as 
Sugimoto-Ikeda-Toki (SIT) interaction, the expectation 
value of the tensor potential amounts to —30 McV and 
the L=2 admixture to 7.3 % in their 4-body state. 

Using the FMD basis with two Gaussians per nucleon, 
in order to allow for a p-wave component in the single- 
particle wave function, we obtain the results presented 
in Tabl e lllli which correspond approximately to those 
in Rcf. [l2|. But the tensor contribution in the FMD 
state is 4.5 MeV larger and the ground state is by 7.25 
MeV more bound, which means that the FMD state is a 
better variational state and represents more tensor corre- 
lations. In contrast to that, as seen in section UlI Al the 
phase-shift equivalent and in this sense realistic potential 
VucOM leads typically to only -5 MeV tensor contribu- 
tion with the same type of FMD trial state. 



TABLE III: Comparison of expectation values of the Hamil- 
tonian, the intrinsic kinetic energy, the total potential energy 
(Coulomb included) and the tensor potential when using the 
trial states of Ref. [I^l, charge and parity projected FMD 
and an exact calculation. For all cases the SIT-potential [12] 
is employed. Numerical values are in MeV. 





Ref. [12] 


FMD-CPP 


Exact 


(H) 


-28.19 


-35.44 


-121.77 


(T) 


64.39 


63.78 


150.56 


(Y) 


-92.58 


-99.22 


-272.33 


(Yt) 


-30.59 


-35.19 


-207.08 



This suggests that the SIT-potcntial has an unrealis- 
tic ratio of tensor to central potential. To investigate 
that further we perform an exact calculation using an 
HH expansion and a Lee-Suzuki transformation to accel- 
erate convergence, as proposed in [43|. The exact results 
shown in Table IIIII exhibit a dramatic overbinding with 
the SIT-potential. The reason is that the tensor potential 
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FIG. 4: (Color online) Off-diagonal momentum space matrix 
elements (in MeV fm^) of the initial AV18 potential (top), 
the VucoM potential for I^ — 0.09 fm^ (center), and the SIT- 
potential (bottom) between the '^5'i and ^Di channel. Dotted 
plane indicates MeV fm^ 



ments of the AVI 8 potential, the corresponding Vucom 
for I^ = 0.09 fm^, and the SIT potential between the ^Si 
and '^Di channel, that is sensitive to the tensor force only, 
are displayed. The off-diagonal matrix elements of the 
initial interaction are strongly reduced for q' ^2 frn" 
by the UCOM transformation Q . In contrast to that the 
SIT matrix elements are not only much larger than those 
of VucoM but even larger than those of the AV18 in- 
teraction and the SIT interaction connects low momenta 
with momenta high above the Fermi momentum. Those 
high momentum components above about 2 fm~^ are 
not present in a low-momentum basis consisting of a few 
Slater determinants and thus the according correlations 
or admixtures cannot be represented. On the other side 
the exact HH state includes high momenta and thus can 
accommodate the correlations caused by this part of the 
potential. This explains the drastic differences seen in 
the energies. 

The SIT-interaction is an example that fitting a phe- 
nomenological potential to a specific nucleus without ref- 
erence to realistic interactions can lead to unstable and 
unpredictable results and does not allow to draw defini- 
tive conclusions on tensor correlations in nuclei. States 
that live in a too small Hilbert space are forced by too 
strong interactions to produce the desired kind of corre- 
lations. Ikeda, Sugimoto and Toki are aware of this prob- 
lem and argue in some of their papers with the strength 
of the G-matrix [10, U^i Il3l- But as Vucom plays the 
same role as a G-matrix we believe that the long range 
tensor correlations can only partly be represented by a 
charge and parity projected intrinsic Slater determinant, 
as shown in Sec. IIII Al 



D. Time reversal symmetry 



In this section we would like to make a remark about 



time reversal symmetry 4J] of our charge and parity pro- 
jected FMD state. With the charge and parity projection 
we actually create intrinsic states which are not invariant 
under time reversal. Since the variational single-particle 
parameters are complex they may get a non-zero imag- 
inary part, violating the time reversal symmetry. In or- 
der to restore the symmetry, we generate a time reversal 
symmetric state by 



*) 



pcMPiixrp 



A/t 



(e'^|Q) + e-'^|Q)) , (33) 



gives a completely unrealistic contribution of -207 MeV. 
Thus neither the FMD nor the trial state of Ref. [i2| can 
represent in a reliable way the huge tensor correlations 
induced by the unrealistic phenomenological tensor part 
of the SIT-interaction. 

The reason for the large difference between the FMD 
or the Sugimoto et al. results |12| and the exact re- 
sult becomes understandable when looking at Fig. 21 
where the non-diagonal momentum space matrix ele- 



where \Q) = T| Q ), is the time reversed Slater determi- 
nant. Then we minimize the energy with respect to the 
phase (j)\ 



Eo 



mm 



(*l#l*) 
(*|*) 



(34) 



We find that for the investigated cases the effect of the 
symmetry restoration is at most of the order of 100 kcV 
and thus negligible for the purpose of this paper. 
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IV. CONCLUSIONS 

Within Fcrmionic Molecular Dynamics the effects of 
mixing proton and neutron components of the single- 
particle states of a single Slater determinant are inves- 
tigated. For that we perform for ^He variational calcu- 
lations by minimizing the energy of the charge and par- 
ity projected Slater determinant using realistic nuclcon- 
nucleon interactions. It turns out that the variation 
needs to be performed after charge and parity projection 
in order to obtain a non-vanishing tensor contribution 
to the ground state energy of the doubly magic nucleus 
"'He. The intrinsic state that is not yet charge projected 
features a non-vanishing pseudo-vector iso-scalar spin- 
isospin density that is intimately related to the pseudo- 
scalar iso- vector pion fields which are responsible for the 
long range part of the tensor force. 

A Hartree-Fock type variation without projection does 
not break charge and parity of the single-particle states 
and hence the expectation value of the tensor interaction 
is zero. 

The extra correlation energy obtained by the new de- 
grees of freedom that mix charge turns out to be small 
for realistic interactions, smaller than anticipated from 



earlier work by Sugimoto, Ikeda and Toki. The main 
reason is that they did not use a realistic interaction but 
adopted a tensor force that scatters to high momenta and 
at low momenta is about twice the strength of the tensor 
part in Vucom- 

Our result is that charge mixing and parity breaking 
of one-body states can account only for a fraction of the 
long range correlation energy missing in a mean-field pic- 
ture (single Slater determinant). One should however 
keep in mind that ^He is the most demanding nucleus in 
this respect. Even if the additional energy due to long 
range correlations is found to be small here, we cannot 
extrapolate this result to open shell nuclei. Also a vari- 
ation after angular momentum projection might change 
the situation. These issues will be subject of future in- 
vestigations. 
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